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ABSTRACT 

As recently shown by Viola et al., the common (KSB) method for measuring weak gravitational shear creates a non-linear relation 
between the measured and the true shear of objects. We investigate here what effect such a non-linear calibration relation may have on 
cosmological parameter estimates from weak lensing if a simpler, linear calibration relation is assumed. We show that the non-linear 
relation introduces a bias in the shear-correlation amplitude and thus a bias in the cosmological parameters Q. m0 and cr 8 . Its direction 
and magnitude depends on whether the point-spread function is narrow or wide compared to the galaxy images from which the shear 
is measured. Substantial over- or underestimates of the cosmological parameters are equally possible, depending also on the variant 
of the KSB method. Our results show that for trustable cosmological-parameter estimates from measurements of weak lensing, one 
must verify that the method employed is free from ellipticity-dependent biases or monitor that the calibration relation inferred from 
simulations is applicable to the survey at hand. 



1. Introduction 

Measurements of weak gravitational lensing have developed into 
one of the main diagnostic tools for the dark-matter distribu- 
tion on large scales. In principle, the ellipticity of the surface- 
brightness distribution in the images of distant galaxies is quan- 
tified in some way, averaged over sufficiently many images, and 
related to the mean ellipticity expected in presence of gravita- 
tional shear to obtain a local estimate of the shear. Even though 
typical lensing effects are very weak and superposed on a sub- 
stantial shape-noise contribution from the galaxies, many mea- 
surements have succeeded and routinely produce cosmological 
parameter estimates which are typically well in agreement with 
alternative determinations, or supplementing them in a highly 
plausible way. See Bartelmann (2010); Bartelmann & Schneider 
(2001) for recent reviews. 

There are two ways of quantifying the light distribution of 
faint galaxies. One of them, a model-free approach, measures 
sufficiently high-order moments of the surface-brightness distri- 
bution and combines them into ellipticity estimates that can then 
be compared to the shear (Kaiser et al., 1995). The other com- 
pares model images to real data and varies the applied shear un- 
til both match optimally (Miller et al., 2007; Bernstein & Jarvis, 
2002; Kuijken, 1999). We are here concerned with the first ap- 
proach, which has the advantage of not assuming an intrinsic im- 
age shape. Voigt & Bridle (2010) have shown how too simplistic 
or rigid model shapes can fundamentally limit the accuracy of 
shear measurements. 

Let 1(0) be the surface-brightness distribution of a galaxy 
image, then the moments are defined as 



Qij...k = J d 2 6I(8)W(6)8 i e j ...6 k , 



(1) 



with the integral is carried out over the image and the weight 
function is introduced to cut the integration off in order to limit 



the inclusion of noise. The three independent second moments 
Qij are combined to form the complex ellipticity 



X : 



gn - 622 + 21612 

Gn + 222 



(2) 



The intrinsic ellipticity of a source,;^, is related to the ellipticity 
X of a sheared image by 



X~2g + g 2 x* 
1 + \g\ 2 - 2%(g X *) ' 



X' = 



(3) 



The reduced shear g = y(l - k)~ 1 appears because the shear y 
itself is measurable only in combination with the convergence k. 
Since y e C, so is g. 

Two essential problems arise in any application of Eq. (2) 
to shear measurements. First, the equation is strictly non-linear. 
The lowest-order linear approximation, 



X s =X~2g, 



(4) 



is typically not sufficiently precise, thus higher-order corrections 
need to be applied. Second, the observed image is further dis- 
torted after being gravitationally sheared. It is convolved with 
the point-spread function of the optical system and of the at- 
mospheric seeing and it is pixellised by the detector. Moreover, 
intrinsically highly elliptical sources are less susceptible to the 
gravitational shear applied. These effects need to be corrected. 

The standard procedure for these corrections has been de- 
fined by Kaiser et al. (1995, hereafter KSB). Although we adopt 
KSB here as a specific example against which we can quantify 
our statements, we emphasise that the central issue of this pa- 
per is by no means a criticism of KSB, but the bias caused by 
any non-linear calibration relation that is approximated by a lin- 
ear one. The KSB method is particularly relevant in this context 
because it is very fast and has been shown to perform well for 
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sources with high signal-to-noise (Bridle et al., 2010). It is rou- 
tinely being used also in shear analyses of galaxy clusters (see 
Hoekstra, 2007, for a recent overview), where calibration biases 
are more severe than for cosmic shear. 

In a recent analysis of KSB (Viola et al., 2010), we have 
identified three potentially problematic steps or assumptions, 
which are: (1) KSB averages shear measurements from galaxy 
images, rather than measuring the shear from averaged galaxy 
images, but the shear measurement and the average do not gen- 
erally commute. (2) KSB implicitly assumes that galaxy ellip- 
ticities are small, while weak gravitational lensing only assures 
that the change in ellipticity due to the shear is small. (3) KSB 
approximately corrects for the convolution with the point-spread 
function (PSF), but does not deconvolve it. Step (1) leads to bi- 
ased results, while assumptions (2) and (3) partially counteract 
in a way dependent on the width of the weight function and of 
the PSF. These effects were analysed in detail in Viola et al. 
(2010). Neither of these assumptions is inherent to the moment- 
based approach to weak lensing. In fact, we recently proposed a 
novel method, dubbed DEIMOS (Melchior et al., 2010), which 
avoids these assumptions and whose performance is therefore 
much more stable against variations of the PSF shape and the 
source ellipticity. 

Here, we study the implications of ellipticity-dependent bi- 
ases, by the specific example of KSB, for the constraints on cos- 
mological parameters. Section 2 summarises the shear calibra- 
tion relations for different variants of KSB. We describe simula- 
tions of the resulting biases in Sect. 3 and present our results in 
Sect. 4. 

2. Shear calibration relations 

The KSB method describes the relation between the ellipticities 
of source and image by a tensor P sh , 



(5) 



Ideally, to lowest order, P is twice the identity, P % ^„ = 25 a p, 
as Eq. (4) shows. The non-linearity of Eq. (3), the presence of a 
weight function W in Eq. ( 1 ) and the convolution by the point- 
spread function complicate matters considerably. Actual imple- 
mentations of the KSB method differ in the approximations and 
assumptions made in the concrete representation of the tensor 
P sh . As in Viola et al. (2010), we investigate four of them here: 

- The original KSB method, labelled KSB, which treats the 
order of the reduced shear and the image ellipticity inconsis- 
tently; 

- A simplification of KSB, labelled KSB-1, which stays con- 
sistently at the linear order; 

- A common modification of KSB, labelled KSB-tr, which ap- 
proximates the tensor P sh by 



P%= Jtr(/> sh )<W 



(6) 



which lacks mathematical justification but performs well in 
practice (see Schrabback et al., 2010, for an application to 
the COSMOS field with its narrow PSF); 
- A consistent extension of KSB to third order in shear and 
ellipticity, labelled KSB-3, which we developed in Viola et al. 
(2010). 

We apply these methods of shear measurement to simulated 
noise-free galaxy images. This proceeds along the following 
steps: 



- A value for the reduced shear g - g\ + igi is randomly drawn 
from a flat distribution for the gi between [0, 0.1]. 

- A value for the intrinsic ellipticity x s is drawn such that 
its modulus [^- s | = 0.3 and its orientation is random within 
[0, 27r). For simplicity, we ignore an intrinsic ellipticity dis- 
persion. 

- Each galaxy is sheared with the reduced shear g. This is re- 
peated with N ~ 100 galaxies to suppress the uncertainty 
due to shape noise in the shear estimate. 

- All galaxies are convolved with a given PSF with a Moffat 
profile. The size of the weight function is set equal to the size 
of the convolved galaxy image. 

- The ellipticity of each image ^ obs is measured. From each 
measurement, a shear estimate g is computed with each of 
the four variants of the KSB method. 



The results are illustrated in Fig. 1 for two different choices of 
the PSF, where the measured reduced-shear estimate g is plotted 
against the true reduced shear g. 
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Fig. 1. The difference between the measured reduced-shear estimate 
and the true reduced shear is shown as a function of the true reduced 
shear, for different measurement methods as indicated by the symbols. 
The ideal case of a measurement obtaining the true value exactly cor- 
responds to the horizontal zero line. Top panel: ground-based obser- 
vations with a wide point-spread function. Bottom panel: space-borne 
observations with a point-spread function of negligible width. 
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Finally, the shear estimates are averaged, accounting for the 
shear responsivity, to arrive at the average shear estimate 



1 



l-X* 2 /2N 



1 N 



(7) 



The results shown in the two panels of Fig.l are based on a 
narrow PSF (bottom) mimicking a space-borne observation, and 
a wide PSF (top) that could be typical for a ground-based ob- 
servation. Both PSFs are fully characterised by their Moffat ex- 
ponent /3 and their FWHM. For the narrow PSF, (J3, FWHM) = 
(2,0.5/? e ); for the wide PSF, OS, FWHM) = (5,5 R e ), where R e 
is the scale radius of the Sersic profile assumed for the galaxy 
images. The ideal relation (following the diagonal) is marked 
by a dash-dotted line. These figures illustrate three main points. 
First, the measured shear estimate falls below or above the true 
shear, depending on the method used for the conversion of ob- 
served ellipticities to shear. Second, the relations between mea- 
sured shear estimate and the true shear are non-linear, which will 
turn into our main issue for this paper. Third, the deviation of the 
shear estimate from the true shear, even its sign, depend on the 
width of the PSF. While KSB-3, for example, performs almost 
perfectly for a narrow PSF, it underestimates the true shear if 
the PSF is wide. The origin of these trends has been identified 
in Viola et al. (2010). Here, we work out the consequences for 
the cosmological-parameter determinations from such measure- 
ments. 



3. Biases on cosmological parameters 

Given the non-linear relations between the shear estimate and 
the true shear, we now consider the following situation. Suppose 
a cosmic-shear observation is being carried out, ellipticities are 
measured and converted to a shear estimate following one of the 
KSB variants described above, the shear correlation function is 
measured and cosmological parameters (essentially the matter 
density £2 m o and the normalisation parameter erg). Suppose fur- 
ther that the shear estimate is in fact calibrated in the analysis 
process, e.g. by means of synthetic data, however assuming a lin- 
ear rather than the underlying, non-linear relationship between 
the shear estimate and the true shear. In other words, the calibra- 
tion is supposed to be carried out assuming that the response of 
the measurement to the shear is independent of the shear itself, 
rather than depending on it. Given such a procedure, which can 
perhaps be considered close to the common practice of identi- 
fying additive and multiplicative shear biases (Heymans et al., 
2006), how biased would the cosmological parameters be, if at 
all? 

In order to answer this question, the following steps need to 
be carried out: 

First, we fix a reference cosmological model and compute 
the shear correlation function for it. To be specific, we choose 
the correlation function g+(6), which is given by 



U0) 



= 2n f 

Jo 



mP K (l)F(lO) , F(x) = 



h(x) 
An 2 



(8) 



The power spectrum P K (J) of the convergence is the usual 
weighted projection of the dark-matter fluctuation power spec- 
trum P$(k) in Limber's approximation. In a spatially-flat uni- 
verse for sources at a fixed redshift, 



PS) 



dw 



w s a(w) 



(9) 



where w and w s are the comoving angular-diameter distances 
to the lensing matter fluctuation and to the sources, respectively 
(see, e.g., Bartelmann 2010; Bartelmann & Schneider 2001 for 
a derivation). 

Second, we choose hypothetical survey parameters, define 
bins {#,} for the measurement of the correlation function and 
compute the covariance matrix C,y for the correlation function 
£ + between these bins. Doing so, we follow the procedure de- 
veloped by Joachimi et al. (2008). This allows the computation 
of two important ingredients. First, we can choose the bins such 
that the signal-to-noise ratio is approximately the same in each 
bin. Second, we can calculate the expected uncertainties of the 
measured correlation function including their mutual correlation 
between bins. The latter step is conveniently achieved by rotat- 
ing the bin vector into the principal-axis frame of the covariance 
matrix, drawing independent Gaussian random numbers with the 
appropriate variance, and inverting the rotation. 

Next, we fit the calibration relations shown in Fig. 1 with a 
realistic quadratic dependence 



g = a + bg + cg 

and an assumed linear dependence 
g'=a' + b'g . 



(10) 



(11) 



These correlation functions of the shear estimate and the true 
shear are then related by 



£ + = a 2 + laccr 2 + b 2 £ + , £ = a' 2 + 



(12) 



and must be corrected accordingly by an offset and an amplitude 
to arrive at £ + . The offsets are 



5%+ - a 2 + laccr 2 and df = a 



a 



(13) 



for the quadratic and linear relations, while the amplitudes are 
b 2 and b' 2 , respectively. We shall assume in the following that 
the offset can be faithfully corrected in any survey by demand- 
ing that the correlation function approach zero at large angles. 
However, the correction of the amplitude remains. Our main is- 
sue for the following discussion is that the amplitude b' 2 derived 
from the assumed linear calibration relation differs from the am- 
plitude b 2 expected from the realistic, quadratic calibration re- 
lation. We thus assume that the measured shear correlation is 
corrected by b'~ 2 while it should be corrected by b~ 2 . That is, we 
erroneously infer the correlation function 



, _ b 2 



(14) 



instead of f+, which is then compared to theory and gives rise to 
a bias in the cosmological parameters. 

Our third and final step is thus quite simple: We take the sim- 
ulated correlation function including its modelled uncertainty 
and error bars, multiply it by the factor b 2 /b' 2 , and then fit theo- 
retical correlation functions to them minimising the x 2 function 



X 



[Cj-€ + (.6i,p)]c^ [Cj-€ + (fij,P)] 



(15) 



where p symbolically abbreviates all parameters entering the 
calculation of the theoretically expected correlation function g+. 

The^f 2 contours following from this procedure are shown for 
three variants of KSB in Fig. 2. 
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Fig. 2. Likelihood contours are shown in the Cl m o-crg plane for: A CFHT-like survey in the top row and a wide survey in the middle and bottom 
rows; three variants of the KSB method, i.e. KSB-tr in the left, KSB-1 in the middle and KSB-3 in the right columns; a wide PSF in the top and 
middle rows and a narrow PSF in the bottom row. The Figure shows that cosmological parameters tend to be biased high for wide and low for 
narrow PSFs, and that the magnitude of the bias depends on the variant of KSB chosen. 



4. Results and conclusion 

The x 1 contours shown in these Figures were calculated as- 
suming two different types of survey, whose parameters con- 
trol the covariance matrices C,y. The assumed survey parame- 
ters are listed in Tab. 1 . The cosmological parameters are set to 
£2m0 = 0.3, Qao = 1 - O m o, h = 0.7 and erg = 0.8 in both cases. 



Name 


CFHT 


wide 


A 


34.2 


100.0 


"gal 


13.3 


30.0 


0V 


0.42 


0.30 


Zs 


0.9 


1.0 



Table 1. Assumed parameters for a CFHT-like and a wider weak- 
lensing survey. A is the area in square degrees, n ga i the number density 
of background galaxies per square arc minute, a € is the variance of the 
intrinsic ellipticity, and z s is the mean source redshift. 



We leave out the original version of KSB because most of 
its likelihood contours fall off the parameter region shown. This 
is because the calibration relation of the original KSB method 
has the largest non-linear contribution and thus produces the 



strongest bias. The main pieces of information evident from 
Fig. 2 are: 

- The wide survey (results shown in the lower two rows) nar- 
rows the contours compared to the CFHT-like survey, but 
does otherwise not affect the results. 

- For a wide PSF (results shown in the upper two rows), cos- 
mological parameters inferred from the likelihood contours 
are biased high for the first- and third-order KSB variants 
(second and third row), while KSB-tr is almost perfectly un- 
biased even with the narrow contours of the wide survey. The 
bias is mild for KSB-1 and substantial for KSB-3. This re- 
flects the shape of the calibration curves for the wide PSF 
as shown in the lower panel of Fig. 1. Since they are convex 
from above, a linear fit to them overestimates the shear and 
thus leads to an overestimate of the shear correlation ampli- 
tude. 

- For a narrow PSF (results shown in the lower row), the in- 
ferred cosmological parameters are biased low for KSB-tr 
and KSB-1, while they are almost unbiased for KSB-3. The 
explanation is similar as above, since now the calibration re- 
lations are convex from below, causing linear fits to underes- 
timate the shear. 

Of course, the amount and the direction of the bias depend on 
the exact parameters chosen, in particular for the size of the PSF 
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Fig. 3. Comparison between the likelihood contours in the fi,„o-cr 8 plane 
for two different types of bias. Left panel: No bias due to shear calibra- 
tion was assumed, but the source-redshift estimate was assumed to be 
off by Az s = 0.05 and Az s = -0.1, as indicated by the line type. The 
right panel is the same as the bottom-left panel in Fig. 2. A wide sur- 
vey was assumed for both panels. The comparison shows that, for a 
wide survey with narrow PSF, the calibration bias can be comparable or 
larger than a realistic bias due to an erroneous source redshift. 

relative to that of the window function, which is in turn adapted 
to the object size. The main conclusions we draw from the results 
in Fig. 2 are instead: 

- The fact that the response of the KSB method to the shear de- 
pends on the shear itself, leading to a non-linear calibration 
relation between the shear estimate and the true shear, gives 
rise to possibly substantial biases in cosmological-parameter 
estimates if a linear calibration relation is assumed. 

- The success of different variants of KSB depends on the 
width of the PSF, and thus on the circumstances of the sur- 
vey. Narrow PSFs tend to cause overestimates which wide 
PSFs tend to cause underestimates. 

- Good results of one variant of KSB, such as KSB-tr for a 
wide PSF, do not allow to conclude that this variant performs 
well under all circumstances. 

The non-linear calibration relation discussed here was not 
detected by the Shear Testing Programme (STEP, Heymans 
et al., 2006) because of its restriction to weak shear. 

Compared to other potential biases in cosmological con- 
straints from weak lensing, the calibration bias can be substan- 
tial. Figure 3 opposes likelihood contours in the Q. m o-cr$ plane 
for a wide survey, in the left panel without calibration bias in the 
shear measurement but a bias in the source redshift, and in the 
right panel without redshift bias but with the calibration bias of 
the KSB-tr method. The comparison shows that the calibration 
bias can be as large as a bias to a source redshift estimate off by 
Az*-0.1. 

Our main concern is not to criticise KSB or its variants, but 
to emphasise that KSB, in all of its variants, causes a non-linear 
relation between the shear estimate and the reduced shear. For 
some variants, the non-linearity is stronger or weaker, depending 
on circumstances, and some variants are mathematically more 
consistent than others. Our main conclusion, instead, fits into one 
statement: Without taking account of its non-linear calibration 
relation, KSB tends to give biased results. A new method avoid- 
ing these problems, called DEIMOS, was proposed by Melchior 
etal. (2010). 

Other methods to estimate the shear may show similar be- 
haviour. Users concerned with highly accurate shear estimation 
should therefore investigate whether the method employed ex- 
hibits an ellipticity-dependent bias, and if so, whether the shape 
and parameters of the inferred calibration relation are applica- 



ble to the survey at hand. This requires that the simulation, from 
which the calibration is to be inferred, mimics the actual survey 
data sufficiently closely, in particular with respect to PSF width 
and the distributions of galaxy size and ellipticity as well as then- 
variations with redshift and signal-to-noise ratio. 
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